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1. Introduction 

Small quantum systems embedded in the continuum of scattering and 
decay channels, are intensely studied in various fields of Physics (nuclear 
physics, atomic and molecular physics, nanoscience, quantum optics, etc.). 
These different open quantum systems (OQS), in spite of their specific fea- 
tures, exhibit certain generic properties such as the segregation of time 
scales [U [21 [3] , the alignment of near-threshold states with decay channels 
[H El E] , the instability of SM eigenstates at the channel threshold [B] , or 
the resonance crossings [3] . All those counterintuitive properties result 
from the non-hermitian nature of an eigenvalue problem in OQSs. 

Resonances are commonly found in various quantum systems. A special 
feature of nuclear OQS problem is the strong configuration mixing which 
has to be handled in both resonant states and non-resonant scattering con- 
tinuum. A suitable framework is provided by the continuum shell model 
(CSM) [3], either the real-energy CSM [101 [31 [11] or the complex-energy 
CSMP21EE3]. 
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One of those salient features of OQSs is related to the eigenvalue degen- 
eracies. In the past, level degeneracies have been much studied in connection 
with avoided crossings in spectra, focusing mainly on the topological struc- 
ture of Hilbert space and associated geometrical phases [HI [15]. Among 
these degeneracies, one finds an exceptional point (EP) [161 El HB] , where 
two Riemann sheets of eigenvalues are entangled by the square-root type of 
singularity. 

In many-body systems, EPs have been studied only in schematic models, 
such as the Lipkin model [19], the interacting boson model [20], and the 3- 
level pairing model [21] from the general class of Richardson- Gaudin models 
[221 I23j . In this paper, we shall present first studies of an EP and its 
consequences in a realistic nuclear CSM. 

2. Shell Model Embedded in the Continuum 

Shell Model Embedded in the Continuum (SMEC) is a recent realization 
of the real-energy CSM |10[lll]. The total function space in this model con- 
sists of two sets: the set of square-integrable functions used in the standard 
nuclear Shell Model (SM) and the set of scattering states into which the SM 
states are embedded. These two sets are defined by solving the Schrodinger 
equation, separately for discrete SM states (the CQS): HsM^i = e\ S <&i 

and for scattering states (the environment): J2c'(E — ^cc'Xe = 0. Here, 
Hsm is the SM Hamiltonian and H cc = Hq + V c is the standard Hamiltonian 
used in coupled-channel (CC) calculations. Channels are determined by the 
motion of unbound particle in a state lj relative to the A — 1 nucleus with 
all nucleons on bounded single-particle (s.p.) orbits in a CQS eigenstate (a 
SM state) ^f' 1 . 

By means of two functions sets : Q = {$f}, V = {(# }j one can define 
the corresponding projection operators Q, P, and the projected Hamilto- 
nians : QHQ = Hqq, PHP = Hpp. Assuming Q + V = Id, one can de- 
termine the third function set by solving the CC equations with the source 
term: u\ +) {E) = G { p\E)H PQ • , where G ( p } (E) = P(E - Hpp^P is 
the Green's function for the motion of a single nucleon in the V subspace, 
E is the total energy of the nucleus A, and HpQ = PHQ. is the 

continuation of a SM wave function <frf into the continuum V . 

Using the three function sets : {$^1, {£e } an d one constructs 

the solution ^> C E = + P^% in the total function space: 

n = Ce + + 4 + \E))(<S>f\(E - H qq (E))-^)(^\Hqp\C c e) (1) 

i,k 
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TLqq(E) stands for the energy-dependent OQS Hamiltonian in Q subspace: 

TCqq(E) = Hqq + HqpG^p\e)Hpq , which includes the coupling to scat- 
tering states and decay channels. In general, this coupling is non-hermitian, 
and depends both on E and Rqp- Eigenfunctions &f of the OQS Hamil- 
tonian TLqq are linear combinations of SM eigenfunctions {3^}. 



2.1. The method of anamenses 

Construction of the many-body basis is a key problem in the SM de- 
scription of OQSs. In Hilbert space, this construction is based on the de- 
termination of a complete s.p. basis which consists of discrete states and 
a scattering continuum. A consistent formulation of nuclear structure and 
reactions becomes possible in CSM/SMEC if s.p. resonances are removed 
from the scattering continuum to be put in the subspace of discrete states 
|24j . This regularization procedure for s.p. resonances is associated with an 
extraction of localized part of s.p. resonances from the scattering continuum 
which leaves only non-resonant scattering states. Recently, this problem has 
been solved by the parameter-free method of anamneses [25J. 

For a given pole of the scattering matrix at k res , scattering states with 

k ~ |Re(/c res ) 2 j in the inner region (r < R) of a basis generating poten- 
tial resemble a bound state wave function. Hence, for a construction of a 
quasi-bound state embedded in a continuum (QBSEC), it is quite natural 
to select the scattering state with a (real) energy e res = H 2 k 2 /2^l. In the 
inner region [0, R], where R is yet arbitrary matching radius for inner and 
outer solutions, the QBSEC should be proportional to the regular solution 
of the Schrodinger equation with k = k. In the external region ([R, oo[), one 
requires that the QBSEC has the bound state asymptotic corresponding to 
k = in. Since the wave number is fixed, both in the inner [0, R] and outer 
[R, oo[ regions, therefore the continuity of the QBSEC wave function and 
its first derivative at r = R yields a unique solution for R, providing an 
unambiguous determination of the QBSEC. 

Extracting such QBSECs from the scattering continuum allows to de- 
termine a new orthogonal subspace of s.p. states, consisting of bound states 
and QBSECs, and a new subspace of non-resonant scattering states. Such a 
s.p. basis is complete and can be directly used for the construction of a com- 
plete many-body basis in Hilbert space. It should be stressed that QBSECs 
constructed in this way yield an anamnesis (a reminder) of resonances in the 
space of £ 2 -functions, removing all resonant aspects from calculated phase 
shifts and providing a reliable continuation of weakly-bound s.p. states 
into a low-energy continuum. This aspect is crucial to achieve a correct 
description of the near-threshold behavior of many-body states. 
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Fig. 1. Eigenenergies of the effective Hamiltonian are plotted for the two lowest 
+ and 2 + eigenstates of 16 Ne. Arrows denote positions of the lowest lp-emission 



thresholds: 
1-15 



15 



15 



F(l/2+)tgip(Zj) (red arrow), F(5/2+) ®p(Zj) (blue arrow) 



F(l/2-) ® p^-; 

text. 



(black arrow). For more details, see the description in the 



2.2. The Hamiltonian 

In this work, we study the one-proton (lp) continuum of 16 Ne to look 
for a possible existence of an EP in a physically relevant region of energies 
and parameters of the SMEC. Details of SMEC calculations can be found 
elsewhere [TUl [TT] . For the effective interaction in Hqq we take the ZBM 
Hamiltonian [26J. The residual couplings between Q and the embedding 
continuum [Hqp) is given by: V\2 = Vq8(ti — V2) and the strength (Vo = 
1100 MeV-fm 3 ) is adjusted to reproduce the experimental width of 5/2^ 
state in 15 F. 

The real parts of eigenvalues of Hqq are shown in Fig. Q] as a function of 
the total energy of 16 Ne. Zero on the energy scale E is fixed at the first lp- 



emission threshold 



15 



J" 



F(l/2 + ) <8> p(Zj) ■ Effective Hamiltonian is energy 



dependent due to the coupling of discrete states of 16 Ne to the embedding 



15 



lp continuum F(i ; C 7r ) ® pb(Zj) which explicitly depends on the total 

energy of 16 Ne. For E > 0, TCqq is complex-symmetric and its eigenvalues 
are complex. Opening of successive lp-emission thresholds (see arrows in 
Fig. [1]) makes a significant change of Of, 2^, 0^ eigenvalues. Coulomb 
interaction shifts systematically the near-threshold minima of eigenvalues 



zakopaneOl printed on November 25, 2008 



■5 



Ei(E) slightly above the corresponding energies of the lp-emission thresh- 
olds. The main contribution to SM eigenstate 0+ comes from a coupling 



to the channel 
1-15 



15 



channel 



F(l/2+) ® p(ls 1/2 ) 
F(5/2+)®p(la 1/2 



In 2f state, the coupling to the 



dominates. In both cases, the dominant 

continuum coupling is in the s-wave. In the excited state, the contri- 
bution of couplings to the ground state (l/2^~) and the first excited state 
(5/2+) is about the same. 



3. The exceptional point in the continuum of Ne 

In contrast to diabolic points, which may arise in hermitian systems 
when two (or more) real parameters are suitably chosen in a Hamiltonian, 
the appearance of EPs is a generic feature of Hamiltonian systems. EP ap- 
pears in the complex g plane of a Hamiltonian H(g) = Ho+gHi, where both 
Hq and Hi are hermitian and [Ho, Hi] ^ 0. The position of degeneracies 
are indicated by the common roots of equations: 

det [H (g) — EI] = ; ^det [H (g) — EI] = (2) 

The maximal number of these roots is M < n(n — 1), where n is the dimen- 
sion of the vector space (the number of many-body states of given quantum 
numbers). Most common degenerate eigenvalues are single-root (EP) and 
double-root (diabolic point or level crossing) solutions of Eqs. ([2]). 

The non-hermiticity of an operator is essential in this context. Whereas 
a hermitian Hamiltonian can be always diagonalized and eigenvalues, even 
if degenerate, always correspond to distinct eigenvectors, in non-hermitian 
Hamiltonian one may find non-trivial Jordan blocks and points in a pa- 
rameter space of the Hamiltonian where both eigenvalues and eigenvectors 
coalesce. These EPs have drastic effects on systems behavior, especially 
concerning geometric phases and adiabatic features. In contrast to DPs, 
which have been intensely studied in nuclear physics [27], the relevance or 
even the existence of EPs has never been demonstrated in a realistic nuclear 
model. Below, we shall present the first study of EPs in the framework of 
a realistic CSM for unbound states of 16 Ne. 

The EPs are searched here in E — Vq plane, where a variation of the 
non-hermitian continuum-coupling term can be fully explored. Fig. [2] shows 
energies and widths of 2g~, l\ resonances as a function of the total energy 
of 16 Ne for a fixed value of the strength Vq = —1357.56 MeV-fm 3 of the 
Hqp coupling term. 2g~, l\ eigenstates and eigenvectors of the effective 
Hamiltonian coalesce at E = 4.765 MeV, forming an EP. At this singularity, 
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Fig. 2. Energy and width of 2j , 2 J eigenstates of an effective Hamiltonian 7Yqq 
in 16 Nc which form an exceptional point at a total energy E — 4.765 MeV. The 
continuum- coupling strength corresponding to this degeneracy is Vq — —1357.56 
MeV-fm 3 . For more details, see the description in the text. 

the dimension D of a vector space decreases (D = n — 1), and eigenvectors 
acquire a non-trivial topological phase if the EP is encircled. 

Contrary to the well-known level repulsion mechanism in hermitian sys- 
tem, the entanglement of eigenfunctions is seen far away from the EP and 
becomes an important source of configuration mixing in the scattering con- 
tinuum. Fig. [3] shows the evolution of real and imaginary parts of the 
quadrupole moment Q2 for five lowest 2f states of 16 Ne. In these calcu- 
lations, we have taken a realistic value of the continuum-coupling strength 
V = -1100 MeV-fm 3 , away from the value (V = -1357.56 MeV-fm 3 ) where 
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Fig. 3. Real and imaginary parts of the quadrupole moment for different 2f eigen- 
states of the effective Hamiltonian in 16 Ne. The calculation has been performed 
for Vo = —1100 MeV-fm 3 . The index V enumerates 2 + states in ascending order 
of energy. 



an EP exists. Imaginary part of Q2 in continuum states, is a measure of 
the uncertainty of a real part. 

One can notice two distinct phenomena: the threshold effect and the 
EP effect. Particle-emission threshold is a branch point in OQSs. In 

Fig. [3j one may notice a strong effect of the coupling to decay channels 
rl5 -1 2+ 

F(5/2 + ) (g> p(/j) for a pair of eigenstates 2^ , 2j • The threshold ef- 
fect is well localized in E and leads to mutually compensating variations 
in (27/|Q 2 |2i") and (2j|Q 2 |2j). The effect of an EP at E = 4.765 MeV 
(V = -1357.56 MeV-fm 3 ) is seen in the eigenstates 2^, 2|. The EP ef- 
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Fig. 4. Real part of the quadrupole moment for 2$ , 2^ eigenstates of an effective 
Hamiltonian TLqq which are involved in the exceptional point at E = 4.765 MeV 
are plotted as a function of E for three values of the continuum-coupling strength: 
V = -1340, -1357.56 (the EP at E = 4.765 MeV), -1400 MeV-fm 3 . For further 
details, see the discussion in the text. 
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feet is weakly localized in energy E and coupling strength Vo, and leads to 
a radical change of 2^, l\ eigenfunctions, as manifested by a complicate 
dependence of (2^|Q2|2^), (2^|<22i2^} on the total energy E. 

The intricate dependence of many-body wave functions in the neighbor- 
hood of an EP can be seen in Fig. U] on the example of ^-dependence of a 
real part of the quadrupole moment ReQ2 in 2%, l\ eigenstates. At a 'crit- 
ical' value of the continuum-coupling strength (Vo = —1357.56 MeV-fm 3 ), 
at which the EP is found, the quadrupole moments (2jj~|Q2|2jj~), ^X\Q,2\^X) 
exhibit a singular behavior, whereas their sum remains finite and smoothly 
changing around the EP. This behavior is a manifestation of a genuine non- 
separability of these two continuum wave functions. It is interesting to no- 
tice, that the EP effect remains even if Vq is changed by 5-10% away from its 
critical value (Vq = —1357.56 MeV-fm 3 ), i.e. the EP effect is rather robust 
and weakly dependent on a fine-tuning of the Hamiltonian parameters. 

4. Summary and conclusions 

In this work, we have demonstrated that the phenomenon of entangle- 
ment of two continuum wave functions by the square-root type of singularity 
can be found for physically relevant parameters of the nuclear many-body 
Hamiltonian at low and moderately high excitation energies. This genuine 
effect in Hamiltonian OQSs is robust and can be seen in a broad range of 
continuum-coupling strengths around the critical value. Moreover, the en- 
tanglement of wave functions remains visible even far away from the exact 
excitation energy of the EP. These features found in realistic SMEC calcu- 
lations allow to hope that the indirect consequences of the EP effect can be 
actually studied in nuclear physics experiments. Further work is needed to 
determine the most useful reaction/structure observables to study manifes- 
tations of the nuclear EP effect. 
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